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■ Abstract 

Dynamical systems, which are equivariant under the action of a non-trivial symme- 
try group, can possess structurally stable heteroclinic cycles. In this paper we study 
stability properties of a class of structurally stable heteroclinic cycles in K™, which we 
call heteroclinic cycles of type Z. It is well-known that a heteroclinic cycle, which is not 
J> ■ asymptotically stable, can attract nevertheless a positive measure set from its neigh- 

bourhood. We say, that an invariant set X is fragmentarily asymptotically stable, if for 
^vq \ any 8 > the measure of its local basin of attraction B$(X) is positive. A local basin of 

attraction Bg(X) is the set of points, such that trajectories starting there remain in the 
00 ■ S- neighbourhood of X for all t > 0, and are attracted by X as t — > oo. Necessary and 

sufficient conditions for fragmentary asymptotic stability are expressed in the terms 
of eigenvalues and eigenvectors of transition matrices. If all transverse eigenvalues of 
linearisations near steady states involved in the cycle are negative, then fragmentary 
asymptotic stability implies asymptotic stability. In the latter case the condition for 
asymptotic stability is that the transition matrices have an eigenvalue larger than one 
in absolute value. Finally, we discuss bifurcations occurring when the conditions for 
asymptotic stability or for fragmentary asymptotic stability are broken. 
Mathematics Subject Classification: 34C37, 37C80, 37G35, 37G40 



1 Introduction 

A smooth dynamical system 

x = /(x), f:W n ^R n (1) 

can possess various kinds of invariant sets - steady states, periodic orbits, tori, heteroclinic 
cycles and strange attractors. Conditions for asymptotic stability and (local) bifurcations 
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of steady states and periodic orbits are well known (see e.g. [8]). For a steady state the 
conditions for stability are formulated in the terms of eigenvalues of the linearisation near 
the steady state. For a periodic orbit they are expressed in the terms of eigenvalues of 
linearisation of the Poincare return map near the periodic orbit. When the conditions for 
stability cease to be satisfied, a bifurcation of the steady state or of the periodic orbit takes 
place. No complete theory for stability and bifurcations of heteroclinic cycles is yet available. 

Let £1, . . . , £ m G M. n be hyperbolic equilibria of flJJ) and Kj : £j —> j = 1, . . . ,m, 
£m+i = £i> De a set of trajectories from £j to The union of the equilibria and the con- 
necting trajectories is called a heteroclinic cycle. Generically heteroclinic cycles are struc- 
turally unstable, because an arbitrary small perturbation of / breaks a connection between 
two saddle steady states. However, the connections can be structurally stable (or robust) 
if the dynamical system has a non-trivial symmetry group and only symmetric perturba- 
tions are considered [H QUI 120], or if the system is constrained to preserve certain invariant 
subspaces [TO] . 

Heteroclinic cycles that are not asymptotically stable can attract a positive measure set 
from its small neighbourhood [21 El El El [131 [T7] . We call such heteroclinic cycles fragmentarily 
asymptotically stable. In earlier papers, several types of stability were employed to describe 
locally attracting, but not asymptotically stable invariant sets: essential asymptotic stability 
(e.a.s.) (different definitions of e.a.s. sets are given by different authors - cf. [TH [12] and 
[3 H]), relative asymptotic stability [U [21], predominant asymptotic stability [16]. If a set is 
stable in any of these senses, then it is fragmentarily asymptotically stable. If a heteroclinic 
cycle is not fragmentarily asymptotically stable, we call it just unstable. 

Asymptotic stability or fragmentary asymptotic stability of structurally stable hetero- 
clinic cycles was considered in a number of papers. A sufficient condition for asymptotic 
stability of heteroclinic cycles is given in [IT]. A heteroclinic cycle is called simple, if all 
eigenvalues of df(£j) are different and the connecting orbits Kj are one-dimensional. Neces- 
sary and sufficient conditions for asymptotic stability of simple homoclinic and heteroclinic 
cycles in R 4 are given in [SI HUB]; necessary and sufficient conditions for fragmentary asymp- 
totic stability of simple heteroclinic cycles in M 4 are given in [TB] (the term fragmentary 
asymptotic stability is not used there). Conditions for asymptotic stability, essential asymp- 
totic stability or relative asymptotic stability for heteroclinic cycles is particular systems are 
presented in [31 HI [5j El LTl El 1121 [171 HE]- In several of these papers, some bifurcations of 
homoclinic and heteroclinic cycles are also studied [51 El El CEO IT%] . 

In the present paper we introduce a class of (structurally stable simple) heteroclinic 
cycles in R n . All simple heteroclinic cycles studied in the papers cited above, except for the 
so-called type A cycles, belong to this class. We call this class type Z heteroclinic cycles. 

For type Z heteroclinic cycles we derive necessary and sufficient conditions for asymp- 
totic stability and fragmentary asymptotic stability. A cycle is fragmentarily asymptotically 
stable, whenever certain inequalities on eigenvalues and eigenvectors of transition matrices 
associated with the cycle are satisfied. If for any j all transverse eigenvalues of df(£j) are 
negative and the cycle is fragmentarily asymptotically stable, then it is asymptotically stable. 
For any of the inequalities determining asymptotic stability or fragmentary asymptotic sta- 
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bility we discuss, a bifurcation of which kind happens, if the inequality ceases to be satisfied 
control parameter is varied. 

2 Definitions 

2.1 Stability 

Denote by $t(x) a trajectory of the system (CQ) starting at point x. For a set X and a number 
e > 0, an e-neighbourhood of X is the set of points satisfying 

B t {X) = {xG R n : d(x, X) < e}. (2) 

Let X be a compact invariant set of ([T]). Denote by Bg(X) its 5-local basin of attraction, 
defined as 

B S (X) = {x G R n : d($ t (x), X) < 5 for any t > and lim d($ t (x), X) = 0}. (3) 

t— >oo 

Definition 1 A compact invariant set X is asymptotically stable, if for any 5 > there 
exists an e > such that 

B e (X) c B 5 (X). 

Definition 2 A compact invariant set X is fragmentarily asymptotically stable, if for any 
5 > there exists a set V, such that V C B$(X) and 

li{V) > 0. 

(Here fi is the Lebesgue measure in MJ 1 .) 

Evidently, if a set is asymptotically stable, then it is fragmentarily asymptotically stable. 
Definition 3 We say, that a setX is unstable , if there exists 5 > 0, such that n(B$(X)) = 0. 

2.2 Heteroclinic cycles 

In this paper we consider dynamical systems that have a group of symmetries, which we 
denote by T. A dynamical system ([1]) is called T-equivariant, where F C O(n), if / : R n — > M. n 
is a T-equivariant vector field, i.e. 

/( 7 x) = 7/(x), for all 7 G T. 

We assume that the group T is finite. 

Let £1, . . . , £ m be hyperbolic equilibria of ([T]) with stable and unstable manifolds W s (£j) 
and W u (£j), respectively, and Kj = W u (£j) fl W s (£j + i) ^ 0, j = 1, . . . , m, £ m+ i = £1, be a 
set of trajectories from £j to £j+i- 
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Definition 4 A heteroclinic cycle is an invariant set X C M. n comprised of a set of equilibria 
{£i> ■ ■ ) im}, and a set of connecting orbits {/c 1; . . . , K m }. 

Recall that for a group T, acting on M n , the isotropy group of the point x G lR n is the 
subgroup 

H x = {7 e T : 72 = x}, 
and a fixed-point subspace of a subgroup £ C T is the linear subspace 

Fix(E) = {x G R n : cxx = x for all a G E}. 

Definition 5 A heteroclinic cycle is structurally stable (or robust), if for any j , 1 < j < m, 
there exists a fixed-point subspace Pj = Fix(Sj), where Sj C T, such that 

• is a sink in Pj; 

We denote = fl P 3 - and the isotropy subgroup of Lj by Tf evidently, ^ G L,-. 
2.3 Eigenspaces, simple cycles and type Z cycles 

For a structurally stable heteroclinic cycle, eigenvalues of df{^j) can be divided into four 
classes HH [TllIS]: 

• Eigenvalues with associated eigenvectors in Lj are called radial. 

• Eigenvalues with associated eigenvectors in Pj-i Lj are called contracting. 

• Eigenvalues with associated eigenvectors in Pj Lj are called expanding. 

• Eigenvalues not belonging to any of the three above classes are called transverse. 

Definition 6 (adapted from fT3f ). We call a robust heteroclinic cycle X G M n \ {0} simple, 
if for any j 

• all eigenvalues of df(^j) are distinct; 

• dim(P j _i Lj) = 1. 

Denote by Pj the orthogonal complement to Pj in IR n . The basis in Pj can be chosen to 
be comprised of contracting and transverse eigenvectors at ^ , or of expanding and transverse 
eigenvectors at 

Definition 7 We call a simple robust heteroclinic cycle X to be of type Z, if for any j 

• dimPj = dimPj+i; 
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• the isotropy subgroup of Pj, Yjj, decomposes P- into 1- dimensional isotypic compo- 
nents. 

Letter Z in the name of the cycle is chosen as an "opposite" one to A: A heteroclinic 
cycle is called to be of type A, if eigenspaces associated with the dominant (having the 
largest real part) contracting eigenvalue of df(£j), the weakest (having the smallest real 
part) transverse eigenvalue of df(£j), the dominant expanding eigenvalue of and the 

weakest transverse eigenvalue of belong to the same £.,- isotypic component [HI HE]. 

This implies that this isotypic component of P- 1 is at least two-dimensional. For type Z 
cycles, by contrast, all isotypic components of Pj- are required to be one-dimensional. 

The condition dim(Pj_i Lj) = 1 implies, that for any j the contracting eigenspace 
at £j is one- dimensional. Together with the condition dimPj = dimPj + i, this implies that 
the dimension of the expanding eigenspace is also one. Denote by n r the number of radial 
eigenvalues and by n t the number of transverse eigenvalues; for a type Z heteroclinic cycle 
n r and n t are the same for all equilibria and n = n r + n t + 2. The radial eigenvalues and the 
associated eigenvectors near £j are denoted by — Tj = — {r^} and vj = {fL}, 1 < I < n r , 
the contracting ones by —Cj and t>j, the expanding ones by Cj and Vj, and the transverse 
ones by tj = {tjj} and = {ujj}, 1 < I < nt, respectively. 

For type Z cycles we can choose a basis in Pj comprised of contracting and transverse 
eigenvectors at £j, or by expanding and transverse eigenvectors at £7+1. Since all isotypic 
components of Pj are one-dimensional and an eigenvector belongs to an isotypic component, 
the basis {t>j +1 , vj +1 } is a permutation of the basis {v^, vj}, possibly accompanied by a 
change of sign of some vectors. Hence, the matrix mapping components of a vector in the 
basis {v?, vj} to components of the vector in the basis {v vj +1 } is a product A±A> , where 
A 7 is a permutation matrix, and A± is a diagonal matrix with elements +1 and —1 on the 
diagonal. 

Following [9l [131 UHL m order to examine stability we construct a Poincare map in the 
vicinity of the cycle. 

2.4 Collection of maps associated with a heteroclinic cycle 

In subsection 12.31 we have given definitions for radial, contracting, expanding and transverse 
eigenvalues of the linearisation df(^j). Let (u, ?),w,z) be local coordinates near £j in the 
basis, where radial eigenvectors come the first (the respective coordinates are u), followed by 
the contracting and the expanding eigenvectors, the transverse eigenvectors being the last. 
Suppose 5 is small. In a 25-neighbourhood of £j, B 2 j-(£j), defined as 

B 2'si.Q = {(u,v,w,z) : max(|u|, |t>|, \w\, |z|) < 25}, 
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system ([T]) can be approximated by the linear syste 

U = — TjU 

V = —CjV 
W = tjW 
Z = tjZ. 

This gives an accurate approximation of the nonlinear system as long as the linear system 
has no low-order resonances. We denote by (u, v,w,z) the scaled coordinates (u, v,w,z) = 
(u,v,w,z)/5. 

Consider a neighbourhood of a steady state £j. Let (u ,t>o) be the point in Pj-i where 
trajectory intersects with the sphere |u| 2 + v 2 — 1, and q be local coordinates in the 
hyperplane tangent to the sphere at the point (u ,fo). Coordinates (u, v) of a point in the 
hyperplane are related to coordinates q as follows: 

(:H"=("H q ' (5) 

where Dj is a n r x (n r + 1) matrix. Some components of uo can vanish, if e.g. Kj-i belongs 
to an invariant subspace in Pj-i- (Note that Pj is not required to be the smallest possible 
subspace.) t> does not vanish, because it is the component in the contracting direction. 

Near ^ we define two crossections of the heteroclinic cycle. One, denoted by Hj in \ 
is an (n — l)-dimensional hyperplane intersecting connection Kj-i at the point (u , t>o, 0, 0); 
coordinates in the hyperplane are (q, w, z). Another one, Hj OUt \ is parallel to the hyperplane 
w = and intersects connection Kj at the point w — 1; coordinates in the hyperplane are 
(u, v,z). Near £j trajectories of system ([1]) can be approximated by a local map (called 
the first return mapj^ (pj : — > H^ ^ relating a point, where a trajectory enters the 
neighbourhood, to the point, where it exits. In the leading order (see dlj and ([S]) ), the local 
map is 

Mi<nhv,M) = ({(^ + ^^s)^ i/e n,^^,{^-'-^}). (6) 

s=l 

The map can be represented as a superposition <pj = Cj 0t V*j ot , where C* ot : R n — > R n ~ l and 
ptot . ftn-i fin rj!^ e ac ^j on Q f ^- ne ma p ptot on ^ ne q coordinates is presented by (jSJ), and 

1 Bclow we assume that all eigenvalues are real. Definition [7] implies that for type Z cycles transverse, 
contracting and expanding eigenvalues are real. Radial eigenvalues can be complex, but this does not change 
our proof significantly. To simplify the proof, we assume that radial eigenvalues are real as well. 

2 If some components tj <s of tj are positive, then the local map is defined for z s satisfying the inequality 
\z s \ < K(l — <5) |k7 | / e J , where if is a constant and S is small (see [9l [T6] ) . However, this restriction is not 
important, because in order to study stability of a cycle we study stability of a fixed point of a collection 
of maps R w — > R , and the maps are defined for all x 6 M. N . Moreover, the local map is defined only for 
particular signs of v and w. To overcome this complication, we consider group orbits of heteroclinic cycles, 
see subsection 12.51 



6 



(4) 



its action on the w and z coordinates is trivial. The map Cj ot is 



/ U N / r j/ej 



/''tot 



v 
w 

V z / 



uw 



to c ^^ I . (7) 



zw 



tj/ej 



Near connection Kj system ([T]) can be approximated by a global map (also called a 

(out) fV(m) 



connecting diffeomorphism) ipj : H^ out ^ — > H 



z i+i 



where superscripts in the notation of components indicate, whether the components are in 
the local basis near £j or near The (n — 1) x (n — 1) matrix 5* ot presents the map 
in local coordinates near £j (i.e., the basis near is the same, as near and the origin is 
shifted to and the matrix v4* ot relates the coordinates in the two local basises). Each 

matrix is comprised of two diagonal blocks (the respective non-diagonal blocks vanish). The 
first n r x n r blocks, and Aj, approximate the maps acting in Pj, and the second blocks, 
Bj and Aj-, the maps acting in Pj~. The matrix Bj is diagonal, because the action of Ej 
decomposes Pj~ into one- dimensional isotypic components by the definition of type Z cycles. 
As discussed in subsection I2.3j Aj- = A^ Aj. 

Denote §j = ipj ° 4>j '■ Hj — > Hj+l, (superpositions of the local and global maps). The 
Poincare map — > for the cycle is the superposition = g m ° ■ ■ ■ ° 9i, for j > 1 
the Poincare maps are constructed similarly: 

9 {]) = 9j-i o . . . o g x o g m o . . . o g jt 

The coordinates (w, z) used to define the maps cjj are independent of q. Hence, we can 
define maps gj, which are restrictions of the maps gj into the (w, z)-subspace: 

g 3 (w,z) = AfA J B,( { ^X } )- 0) 

We call the set of maps {g™} = {gi, ■ ■ ■ ,g m }, where g<j : M rat+1 — > M rat+1 have been con- 
structed above, a collection of maps, associated with the heteroclinic cycle . . . , £ m }. The 
collection of maps {g 1 ^ 1 } = {gi, • • • , g m , gi, ■ ■ ■ , gi-i} is associated with the heteroclinic cycle 
{£,1, ■ ■ ■ , Cm; 6i> ■ ■ ■ > 6-i}> which geometrically coincides with the former cycle. 

2.5 Comments on definitions of a heteroclinic cycle 

Note that a heteroclinic cycle, as defined by the commonly used definition HJ is never asymp- 
totically stable according to definition [1] due to the following reasons. Denote by 7 (out ) a 
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symmetry satisfying ^ out ^> e Tj and 7 (out ) ^ The symmetry fixes reverses the sign 

of w, maps i7- into itself, and maps the connection kj into a different connection 
exiting a neighbourhood of £j via ry{out) jji° ut ) (which differs from Hj OUt ^ by the sign of u>). 
Hence, for any small e we can find points in B e (X) fl such that trajectories starting 

there do not follow the connection Kj. Similarly, the sign of v in Hj is either negative 
or positive, and the symmetry 7^ G Tj, such that 7^ ^ £j, maps Hj into a different 
crossection with an opposite sign of v and the connection Kj-\ into Hence, the 

map <f)j : Hj — > Hj*' is defined for particular signs of w and v. 

For this reason, a heteroclinic cycle is often defined as (a connected component of) an 
orbit under the action of the group of symmetries T of a heteroclinic cycle in the sense of 
definition HI For a group orbit the maps 4>j are defined for all signs of w and v , because in 
addition to Kj_\ and Kj the orbit involves ^ %r ^Kj_\ and ^ out ^Kj as well. A group orbit can be 
asymptotically stable according to definition [TJ Below, when speaking about a heteroclinic 
cycle, we always assume a group orbit. 

Consider a sample trajectory close to a heteroclinic cycle. After the trajectory passes 
near it can head for following connection Kj, or for ^ out ^j + i following connection 
^(out) K ,^ choice depends on the sign of w. Hence, different trajectories can visit different 
sets of equilibria. Suppose that all subspaces Pj are maximal, i.e. there does not exist an 
invariant proper subspace of Pj, Pj, such that Kj C Pj. Then generically none of the (u ,fo) 
components vanish. The signs of components of z are preserved by the local maps. Any 
global map preserves the sign of each components of z for all trajectories, or reverses it for 
all trajectories. Hence, a particular path along the cycle followed by the trajectory (i.e. the 
sequence of equilibria visited by the trajectory) $t(x) for x e is uniquely determined 

by the signs of coordinates of x (i.e. all points in each orthant of M n_1 of Hj follow the 
same path.) 

Often when definition H]of a heteroclinic cycle is used, a different object is tacitly assumed. 
Suppose there exists a symmetry 7 e T, such that the heteroclinic cycle {£1, . . . , £ m } can be 
generated from its smaller subcycle {£1, . . . , £/}, m = IK, by applying symmetry 7: 

7£ s+(fc _i)z = £ s+k i for all 1 < s < I and 1 < k < K. 

Then sometimes the subcycle {£1, ...,£/} is called a heteroclinic cycle referring to its group 
orbit. If I — 1, then the heteroclinic cycle is called a homoclinic cycle and the symmetry 7 
is called a twist [21 [20] . 

2.6 Collection of maps: definitions of stability 

For a collection of maps {g™} = {gi, . . . ,g m } we define superpositions 

g {]) = g 3 -l o . . . o gi o g m o . . . o g j+1 o 3j (10) 
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and 

{gi°---°9j, l>3 
(11) 
gio...og m og 1 o...og j , l<j 

Given a collection of maps {g™} = {gi, . . . ,g m }, gj : M. N — > M. N , we can define a discrete 
dynamical system as follows: y n +i = g^'Yn- We call y 1 G M. N a fixed point of the collection 
of maps {g™}, if g^y 1 = y 1 - Evidently, y' = gy-i^y 1 is then a fixed point of the collection 
of maps {gl^ 1 }- 

Definition 8 We say, that a fixed point y 1 G M. N of a collection of maps 

{g?} = {gi,---,9m}, g r .R N ^R N , (12) 

is asymptotically stable, if for any 5 > there exists an e > 0, such that for any 1 < I < m 

d(x,y l ) < e, where y l = gy-i^y 1 , implies 
d((g ij) ) k g(j-i,i)X,gu-i,i)y l ) <S for all l<j<m, k > 

and 

lim d((5f 0) ) fe 5f (i _ij ) x,5( (i _ 1 . y') = for all 1 < j < m. 



Definition 9 We say that a fixed pointy 1 G K o/a collection of maps {g™} is fragmentarily 
asymptotically stable, if for any 5 > t/iere exists V C ^({g] 11 }, y 1 ), where 

Bsdg?}^ 1 ) := {x : x G K", rf((g«)^_ ljl)X , ^-i^y 1 ) < 5 for all 1 < I < m, k >0 

and lim d{(g^) k gn- lt i)X, ga-x^y 1 ) = for all 1 < I < m}, 

sttc/i i/iai /^(V) > 0. 

Definition 10 We say that a fixed point y 1 G K.^ of a collection of maps {g}™ is unstable , 
if there exists 5 > 0, such that 

M^({^},y 1 )) = o. 



3 Stability of a cycle and a collection of maps 

In this section we prove two theorems stating that a type Z heteroclinic cycle is (fragmentar- 
ily) asymptotically stable if and only if the fixed point (w, z) = of the collection of maps, 
associated with the cycle, is (fragmentarily) asymptotically stable. (The point (w, z) = is 
evidently a fixed point of the collection of maps constructed in subsection 12.41 ) 
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Theorem 1 Let {g™}, gj '■ M n ' +1 — > M nf+1 , be a collection of maps associated with a hete- 
roclinic cycle of type Z. Then the cycle is asymptotically stable, if and only if the fixed point 
(w, z) = of the collection of maps is asymptotically stable. 

Proof: A necessary condition for asymptotic stability of both the cycle and the collection of 
maps is that all transverse eigenvalues are negative. From now on, in this proof we assume 
that this is the case. 

For a trajectory $t(x) belonging to B$(X) for all t > 0, denote by $^(x) the A;-th 
intersection of $t(x) with Hj , by _^-°^(x) the k-th intersection of $t(x) with Hj OUt \ and 

by tjk an d tjk the times, when the intersections take place. 

For a sufficiently small 5 > the collection of maps {g™}, where gj : M n_1 — > IR."" 1 , 
accurately approximates trajectories in the vicinity of the cycle; i.e., if x G H^ and $t(x) G 
B 5 (X) for all t > 0, then 

« is^h-x^ (13) 

Stability of the cycle implies stability of the fixed point of the collection of maps, because 

• if $t(x) G B$(X) for any x satisfying „(x, X) < e, then 

d($$(x),X) < 5 for any j, k and x G Hf n \ |x| < e, 

and hence, by (HH, \(g^) k gj-i^\ < \(g^) k g^ hl x\ < 5; 

• if Hindoo <i(_ , t (x), X) = for any x satisfying „(x, X) < e, then 



lim dml\x),X) = for any j and x G H\ m \ |x| < e, 



->oo 

and hence, by (fT3|) . 



lim |(j«)Vi,i x I < , lim |(_ 0) )%-mx| =0. 

fe— >oo /c — ^oo 

Therefore, by definition [HI is an asymptotically stable fixed point of the collection of maps. 

To prove that asymptotic stability of the cycle follows from asymptotic stability of the 
fixed point of the collection of maps, we first show, that asymptotic stability of the 
fixed point (w, z) = of the collection {g™} implies asymptotic stability of the fixed point 
(q, w,z) = of the collection {g™}- Denote by gj the map gj restricted to the subspace 
(w — 0, z = 0) equipped with the q coordinates. By virtue of j_])-(j_J and ([8]), the map takes 
the form 

g q j (q,w,z)=A»B»C»Vf t [ I |, (14) 
/ u \ 



i } V _) ~j -y ^3 3 3 J 

where 



v 
w 

V z / 



({ Ul w r ^}). (15) 
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Let K j A = n r m.axi }S \Aj >la \, K j B = n r maxi )S \Bj tla \ K J D = n r max(|u |, 1, max /jS |-D|j) and 
r J min = mxD.ifrjj/ej). For a given 5 > 0, choose a Sj satisfying 

K^Kj.Kpd'j' < 5/2. (16) 
Since is a stable fixed point of the collection {g™}, we can find an e > such that 

\(gW) k gj_uX.\ < min(5/2,min(5 s ) for all 1 < j, I < m, k > 0, |x| < e. (17) 

s 

If ( fT7l) holds true, then 

\{g Q) ) k 9j-i,iA < 8 for all 1 < j, I < m, k>0, |x| < e (18) 

by virtue of ©-(E]) and ffl4j) - ffl6]) . The proof, that lim fc ^ oo (^ ( ' J ' ) ) fc ^ (i _i i0 x = implies 
lim / ! C _ s>00 (^^^) fc ^ :; _i i ;x = 0, is similar and it is omitted here. 

Second, we prove that if the fixed point of the collection {g™} is asymptotically stable, 
then for any 5 > we can find an e > 0, such that for any x £ Hi and |x| < e the 
trajectories $t(x) satisfy 

d($t(x),X) < 5 for any t > 0. 

In a sufficiently small neighbourhood of the origin at the hyperplane H^ \ the map 
ipj : flj^f — > Hj is predominantly linear. For any intermediate crossection H^ between 
Btf and BT the induced ,nap fff <> -> *f is a ls o predo— y linea, Hence, we 
can find a positive _ft"j glob - ) such that 

d($ t (x), X) < Kf oh) d(^;\x),X) for any f, < t < t<$> . 

Now consider the trajectory $t(x) at the time interval < t < Near we 

project the cycle X and the trajectory $t(x) onto the plane (v,w) and onto the orthogonal 
hyperplane (u, z). By d^"' w \-, •) and S u ' z ^(-, •) we denote the distances between the projec- 
tions onto the (v, w) plane and onto the (u, z) hyperplane, respectively. Simple algebra (not 
presented here) attests that 

d("' w >($ t (x),x) < (d^f;\^),x)p/^\ 

where {3j = Cj/ej. The estimate 

^ z )($,(x),x) < d(*$(x),x) + d($J.;r ) (x),x). 

follows from (j3J). We denote 

^(giob) = max (^(g lob )) ; ^ = max (fl.). 

A J A 
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For a given S > 0, suppose Si > satisfies 

K (g\oh) 5 P/(l+P) < 5 / 2; (^(glob) + ^ < ^ 

Since is asymptotically stable, we can find e > such that d(Q^\x), X) < Si for all j 

and k, provided x G for some I and d(x, X) < e. The estimates presented above imply 
that d($ t (x),X) < S for all t > 0. 

Finally, we consider x H^ for all /. Denote by x the first intersection of $t(x) with 
some Hi . By the arguments similar to those presented above, at least one of the estimates 

d(5t,X) < Kd(pc,X), if x is near k\ 

and 

d(x, X) < X(d(x, X) + d(x, Xf /{l+l3) ), if x is near 

holds true. Each of the two inequalities imply, that for e defined in the previous paragraph 
we can find an e > 0, such that d(x, X) < e implies d(x, X) < e. Hence, for d(x, X) < e the 
estimate d($ t (x),X) < 5 holds true for all t > 0. 

The proof that lim^oo d($ t (x),X) = follows from limfc_ 5>00 (^^^) fc ^(j_i i ;)X = is similar 
and is omitted here. QED 



Lemma 1 Let X be a type Z heteroclinic cycle. If X is fragmentarily asymptotically stable, 
then for any S > and any j 

/1 "- 1 (ff ) n8 i (i))>o, (19) 

where /i n_1 is the Lebesgue measure in IR n_1 . 

Proof: Denote Qj = H^ PI ^(X). Suppose (fT9l) is not satisfied for some j, i.e. 

/ i"- 1 (g,) = o. (20) 

The set Bs (X) can be regarded as the union of segments of trajectories going from Qj in the 
direction of negative t till an intersection either with ifj m ^ or with the boundary of B$(X) 
occurs: 

B s = {$ t (x) : x G Qj, t final < t < where $ t a na i(x) G or d($ tflna i (x) , X) = S}. 

Denote by £j;f 6) (X) the part ofB s {X) bounded by Hf ut) and Hjj% and by #g c) (X) the part 

bounded by Hj and H^^f'. By linearity of the global map, (J20D implies /i(£j 9 _^j(X)) = 
and 

^(hI? nB 5 (I)) = 0. (21) 
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Since the local map satisfies (j4]) and due to (12T|) . 

MBftW) = and /i"" 1 ^ n B S (X)) = 0. 

Applying the same arguments to the sets Qj-i = H^™\ fl B$(X), Qj-2, ■ m — 1 times, 
we obtain that fi(Bs{X)) = 0, in contradiction with the statement of the lemma. Therefore, 
the assumption ([20]) is false, i.e. ^ n ~ l (Qj) > for all j. QED 

Theorem 2 Let {g™}, gj '■ lR nt+1 — > R nt+1 ; be the collection of maps associated with a type 
Z heteroclinic cycle. Then the cycle is fragmentarily asymptotically stable, if and only if 
is a fragmentarily asymptotically stable fixed point of the collection of maps {g™}- 

Proof: By Lemma 1, for any j ffT9l) holds true. Therefore, the measure /i" t+1 (in R" t+1 ) of 
the orthogonal projection of the set ifj m ^ C\Bs{X) into the plane q = is positive. For a small 
S the collection of maps gives accurate predictions for trajectories $*(x), $t(x) C B$(X) for 
t > 0, and the coordinates w and z are independent of q. Hence, is a fragmentarily 
asymptotically stable fixed point of the collection of maps. 

The proof, that fragmentary asymptotic stability of of the collection of maps {g™} 
implies fragmentary asymptotic stability of of the collection of maps {g™}, is similar to 
the one presented in the proof of theorem [1] and is omitted here. 

Let the constants K( gloh \ K and (3 be defined as in the proof of theorem [TJ The same 
arguments as employed in this proof imply that the inequalities 

d($ t (x),X) < K^d^\ic),X) for any j, k and t, t£$ k < t < t$ if d($£°(x), X) < 5; 

d^($ t (x),X) < (d($%\x),X))M 1+ ® for any j, k and t, t<$> < t < iff if d{*<$\x),X) < 5; 
^ z )($,(x),X) < d(*${x),X) + d{*«f{x),X) for any j, k and t, t$ < t < tff\ 

if d(*%\x),X) < 5 and d($£f >(x), X) < 5. 

are satisfied for a sufficiently small 5 > 0. 

Defining 5± as in the proof of theorem [H we find that <i($ t (x),X) < 5 for all t > 0, 
provided x G Hi is such that x G ^({^x 1 }, 0). The proof, that the condition x G 
Bs^ig™}, 0) implies Hindoo d($ t (x), X) = 0, is similar. 

Let h[ 1 ^ and _ a be two hyperplanes, parallel to Hi and located at distance s from 
this hyperplane. Consider the set Qi comprised of pieces of trajectories, contained between 
the hyperplanes H^ and H[ zr ^ s , whose points of intersection with the hyperplane Hi 
constitute the set Qi := B Sl ({g™},0). For small s, Qi C Bs{X) and fi(Qi) = 2sfj, n ~ 1 (Q j ) > 
0. Thus, the set X is fragmentarily asymptotically stable. QED 
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4 Stability of fixed points of a collection of maps 



4.1 Transition matrix 

Denote by Aij the maps gj in the new coordinates^) rj, where 

rj = (ln|iw|,ln|zi|, ...,ln|z IH |). (22) 
As discussed in subsection I2.4[ the maps are linear and have the structure 

Aijtj = Mjtj + Fj, (23) 

where 





( hi 





. 


• °\ 




ha 


1 


. 


. 


M S := AjBj = Aj 


hj 





1 . 


. 









. 





are basic transition matrices of the maps. Here Aj and Bj are N x N matrices, N = nt + 1, 
Aj is a permutation matrix and the entries bjj of the matrix Bj depend on the eigenvalues 
of the linearisation df(£j) of (p]) near £j as follows: 

bj t i = Cjjtj and &j,z+i = —tj,i/ej, 1 < I < nt, 1 < j < m. (25) 

We call {Ai™}, like {g™}, a collection of maps, associated with the heteroclinic cycle. A 
fixed point (w, z) = of the collection {g™} becomes a fixed point r) = — oo of the collection 
{A4™}. In the study of stability of the point (w, z) = we consider asymptotically small z 
and w, i.e., asymptotically large negative 77, and hence finite Fj can be ignored. 

Transition matrices of the superposition of maps g^' and gjj are the products M^' = 
Mj^ x ■ ■ ■ M\M m ■ ■ ■ M j+1 Mj and M jtl = Mj ■ ■ ■ M x M m ■ ■ ■ M i+1 Mi (or Mj ■ ■ ■ M t if j > I), 
respectively. For a collection of permutation matrices Aj we define similarly Aw and Ajj. 
Denote by A s eigenvalues of matrices M w (they are independent of j, since all matrices Mw 
are similar) assuming the descending ordering of their real parts (generically, the real parts 
are all distinct, except for pairs of complex conjugate eigenvalues). Let also w^' s denote the 
eigenvector of the matrix associated with the eigenvalue \ s . 

4.2 Two types of eigenvalues of a transition matrix 

Consider a matrix M := = M m . . .M\ : M. N — > 1H N , which is a product of basic 

transition matrices of the form ff24l) . We separate the coordinate vectors ej, 1 < / < N, 
into two groups. The first group is comprised of the vectors e/, for which there exist such 

3 Here we ignore the matrix A±, because it is irrelevant in the study of stability. It becomes important 
in the study of bifurcations - A± determines the length of periodic orbit (s) bifurcating from a heteroclinic 
cycle (see subsection 15. 2[) . 
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k and j that (A^) k Aj-i^ei = ei (recall that Aj are permutation matrices), the second one 
incorporates the remaining vectors. Denote by V slg the subspace spanned by vectors from 
the first group, and by V ms - the subspace spanned by vectors from the second group (the 
superscript "ins" and "sig" stand for significant and insignificant, respectively). 

Theorem 3 Let V S1S and V ms be the subspaces defined above. 

(a) The subspace V ms is M-invariant and the absolute value of all eigenvalues associated 
with eigenvectors from this subspace is one; 

(b) Generically all components of eigenvectors, not belonging to V ms , are non-zero. 

Proof: (a) Denote by V^ ms , 2 < j < m, the subspaces constructed for matrices similarly 
to V ins . Since e x £ 1A ins , the action of Bj on Vj™ is trivial. Evidently, AjVj m = Vjl\. Hence, 
the actions of M and A coincide on the subspace V ms . A-invariance of V ms follows from the 
definition of V ms , and hence V ms is M-invariant. The matrix A is a permutation matrix, 
because it is a product of permutation matrices, and thus all its eigenvalues have the unit 
absolute value. Therefore, this holds true for the restriction of M on V ms . 

(b) Let w 1 ' 5 be an eigenvector of M, which does not belong to V ms . At least one significant 
component of w 1,s does not vanish. Denote it by wj' s . By definition of V S1S , there exist such 
j and k that (A^) k Aj-i t xei = e%. Eigenvectors w J ' s and w 1 ' 5 are related as follows: 

w j,s = (M^) fc M i _i >1 w 1 ' a 

(up to an arbitrary factor). Hence, generically the first component of w J,s does not vanish. 
Since eigenvectors satisfy the relation 

w - = M m>j w j > s 

and holds true for all matrices My in the product M mj , generically all components of 
w 1 ' 5 are non-zero. QED 

We call insignificant the eigenvalues associated with eigenvectors from V ms , and signifi- 
cant the rest ones. Generically the absolute values of all significant eigenvalues differ from 
one. 

4.3 Two lemmas 

We consider here a criterion for positivity of the measure of the set 

U-°°(M) = {y: y G R N , lim M k y = -oo} 

k— >OD 

in the terms of the dominant eigenvalue and the associated eigenvector of matrix M : R — > 
R . We have denoted 

K_ = {y = (yi, • • • ,Vn) ■ Vj < for all j}, R N = {y = {y u ...,y N ): Vj < for all j}, 
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U s = {y : max?/, < S}, U s = {y : max?/, < £}. 

s s 

After the change of variables (1221) the Lebesgue measure of a set of initially a finite 
measure can become infinite. To avoid this, the measure of a set V is understood as its 
measure in the original variables. Note that n(V) is strictly positive, if and only if this is 
true for the image of V under the mapping y — > e y , which is the inverse of (12"21) . 

We denote by A max arid w max the maximal in absolute value eigenvalue of the matrix M : 
M. K and the associated eigenvector, respectively. If |A max | > 1, all other eigenvalues 
are supposed to be strictly smaller in absolute value than A r 



>max- 



Lemma 2 fi(U °°(M)) depends on A max and w max as follows: 
(i) //|A max | < I, then C/-°°(M) = 0. 
(it) //A max is real and A max < — 1, then fjL(U~°°(M)) = 0. 
(Hi) If X m3X is complex and |A max | > 1, then n(U~°°(M)) = 0. 

(iv) If A max real, A max > 1 and yj ; max w™ x < /or some I and q, then fi(U~°°(M)) = 0. 

(v) i/A max ?s real, A max > 1, ^ max M;™ ax > for alll<l,q< N, then fi(U~°°(M)) > 0. 
Proof: Consider expansion of y G 

A 1 



^ a-iWi (26) 



i=l 

in the basis comprised of eigenvectors of M, Wj, 1 < z < iV. The k-th iterate is 

N 



M k y = J2^^- (27) 



i=i 



Consequently, 

(1) If |A max | < 1, then 



lim M — oo for any y G 



because if limfc^oo Af exists, it is zero or finite. 
(2) If A max is real and A max < — 1, then 

for any y G R such that a max 7^ the limit lim M k y does not exist, 



k— >oo 



because for k — > 00 the iterates M h y become aligned with A max w max and the signs of 
individual components of A max w max alternate for odd and even k. (Here a max denotes the 
coefficient in front of w max in the sum (1261).') 
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(3) If A max is complex and |A max | > 1, then 

for any y G M. N such that ai imax 7^ or a2,max 7^ 0, the limit lim M k y does not exist, 

because for k — > 00 the iterates M k y are attracted by the plane spanned by w 1,niax and 
w 2 ' max . Only a quadrant of the plane belongs to Mr. The map M expressed in polar 
coordinates in the plane is a multiplication by re 1 ^ (i.e. it involves rotation by angle ip with 
ip 7^ 2ir and multiplication by r > 1). Hence, the limit does not exist. 

(4) If A max is real, A max > 1 and t/j^io" 1331 < for some I and q, then 

for any y such that a max 7^ lim M k y ^ R^. 

k— >oo 

To show this note that for k — > 00 the iterates M k y become asymptotically close to 
a max A^ ax w max , the signs of two individual components, a^X^wf 1 ^ and a^X^w™ 3 *, 
are opposite and hence one of them is positive (unless one or both of the two components 
vanish) . 

(5) If A max is real and A max > 1, ^ max ^m ax > for all I and q (for the sake of definiteness we 
assume W7™ ax > 0), then 

lim M n y = —00 for any y G with a max < 0, 

n— >oo 

because for k — > 00 the iterates M k y become asymptotically close to a max A^ iax w max . 

Evidently, (1), (2), (3) and (4) imply (i), (ii), (hi) and (iv), respectively. 

To prove (v), note that the point — w max belongs to 1R^. Therefore there exists a neigh- 
bourhood V of — w max , such that V C IR^, and for y £ V the sign of a max in the expansion 
(|26|) is negative. The measure of this neighbourhood is positive. Thus (5) implies (v). QED 



Lemma 3 Let M be a matrix with non-negative entries, |A max | > 1 and w ; max 7^ for all 
1 < I < N. Then 

(i) A max is real and positive. 

(11) ffi^w™ > for all I <l,q < N. 

(m) U-°°(M) = R N . 

Proof: Consider an y G Mr. Since all entries of M are non-negative, 

M k (y) G for all k. (28) 

Consider expansion (|26|) for y. If A max is complex or real negative, then some iterates M k y 
are not in ]R_ , as noted in the proof of Lemma [2j Similarly, if wf^^w 1 ? 3 * < for some / and q, 
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then ([2811 does not hold true for sufficiently large k. Since no components of the eigenvector 

w max vari i s h ; ^max^max _/_ g £ Qr a JJ ^ an( J ^ an( J |; nere f ore (jj) jg satisfied. 

To prove (iii), note that if a max > in the expansion ( )26|) for y, then lim^oo M k y = oo, 
which is prohibited by ( 128|) . Next, we show that a max can not vanish. Suppose a max = 
in expansion (1261) for some y G K._. There exists a neig hbourhood C/ C R N of y. We can 
find y in this neighbourhood such that in expansion (1261) for y the factor in front of w max is 
positive, which contradicts fl28|) . Hence, a max < in expansion (T26D for all yeR?. Finally, 
a max < implies \iva k ^ t00 M k y = — oo. QED 



4.4 Properties of maps 

In this subsection we prove two theorems concerning linear maps, which we use to study 
stability of fixed points of a collection of maps. The following lemma is a restatement of 
Lemma [2J 

Lemma 4 Suppose A max is the largest in absolute value eigenvalue of matrix M and w max 
is the associated eigenvector. fi(U~°°{M)) > 0, if and only if the three following conditions 
are satisfied: 

(i) A max is real; 

(H) A max > I; 

(iii) tu™ ax M;™ ax > for all I and q, 1 < I, q < N. 

Evidently, if n(U~°°(M^)) = for some j, where is the transition matrix of the 
collection of maps {gj -1 }, then there exists a small 5 > such that fi(Bs({g™}, 0)) = 0, and 
therefore is an unstable fixed point of the collection {g™}- 

Theorem 4 Let Mj be basic transition matrices of a collection of maps {g™} associated with 
a heteroclinic cycle of type Z. Suppose that for all j, 1 < j < m, all transverse eigenvalues 
of dfi^j) are negative. Then: 

(a) If the inequality |A max | > 1 holds true for the transition matrix M := = M m . . . M 1; 
then is an asymptotically stable fixed point of the collection of maps {g™}. 

(b) If |A max | < 1, then is unstable. 

Proof: (a) The matrices M^') are similar, hence if the maximum absolute value of eigenval- 
ues of matrix M is larger than unity, this is also the case for M^) for any j. All transverse 
eigenvalues of df{^j) being negative, (124"I) and (125!) imply that all entries of matrices Mj are 
non-negative. Since |A max | > 1, Theorem [3] implies that w™ ax 7^ for all 1 < q < N . Hence, 
by Lemma [3] 

lim {M^fM^iy = -00 (29) 

fc— >oo 
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for all 1 < j, I < m and any y G Mr. 

For some S < denote by V ± (S) the intersection of the set Ug with the N—l- dimensional 
hyperplane, orthogonal to the line (1, . . . , 1) and crossing the line at the point (2S, . . . , 2S). 
All entries of matrices Mj are non-negative and the set V A_L (S') is compact, and hence (1291) 
implies that for any / and j there exists a constant R l 'i < such that 

{M^fMj^iy < R lJ (30) 

for all k > and y G V^i^S). Denote 

-Rmax = max R l ' J . 

l<l,j<N 

For a given R < set S = 2SR/ R majX . Linearity of the maps Mj and fl30l implies 

(M^) & M,_ lj/y < i? 

for all k > 0, 1 < j, I < m and y G C/5. 

Changing coordinates from 77 back to (w,z), we conclude that by definition [8] the point 
is an asymptotically stable fixed point of the collection {g™}. 

(b) The statement of the theorem follows from Lemma HI QED 



Theorem 5 Let Mj be basic transition matrices of a collection of maps {g™} associated 
with a heteroclinic cycle of type Z. (For type Z heteroclinic cycles the matrices are of the 
form (24\)-) Denote by j = ji, . . .ji the indices, for which some entries of Mj are negative; 
the entries are all non-negative for all remaining j . Assume L > (the case L = is treated 
by Theorem^. 

(a) If for at least one j — ji + 1 the matrix M^> does not satisfy conditions (i)-(iii) of 
Lemma^ then is an unstable fixed point of the collection {g™}- 

(b) If the matrices M^' satisfy conditions (i)-(iii) of Lemma\4\for all j =ji + l, then is 
fragmentarily asymptotically stable fixed point of the collection {g™}- 

Proof: (a) Matrices M^' satisfy, or not satisfy, conditions (i)-(ii) of Lemma H] simultaneously 
for all j, because all M^' are similar. Hence, by Lemma HI n(U~°°(M^)) = for all j. 
Suppose condition (iii) is not satisfied for some j = J. The iterates (M^) k Mj 

— i.iy become 

aligned with w max,J in the limit k — > 00 (see fl27|) ). Since w™ ax,J w™ ax ' J < for some / and 
q, the iterates escape from Us for any S < 0, when k is sufficiently large. Hence part (a) is 
proved. 

(b) Suppose for some / the matrix M^ l+1 ^ satisfies condition (iii) of Lemma |U Matrices 
Mj, ji + 1 < j < ji + i — 1, have non- negative entries, they are of the form f)24p and w maxJ = 
Mj-i . . . M jl+2 M jl+1 w max ^ +1 . Thus if w ^t+ 1 w ***ji+i > for all I and q, 1 < l,q < m, 
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then w i max '-'w™ ax ' : ' > for all I, q and j, ji + 2 < j < Hence, it suffices to check condition 
(iii) only for j = ji + 1, 1 < I < L. 

Denote by V compl the N — 1-dimensional M- invariant complement to w max in M. N , and 
by Sf _1 (0) the unit ball in V compl centred at 0. If y G Sf _1 (0), then 



lhn ( -^—M^ y - 0. 



Therefore, there exists a constant K\ such that 

1 \ k 

-Ml y < K\ 



.A, 

for all y G i?j v_1 (0) and k > 0. Similarly, there exist constants Kj, j = 2, . . . , JV, such that 

\ '"max / 

for all y G Sf -1 (0) and k > 0. Denote K max = rnaxj ii'j. Let d m i n be the minimal distance 
from the endpoints of eigenvectors w maxJ starting at the origin (we assume now that upon 
normalisation w ma ^ are unit eigenvectors) to iV — 1-dimensional coordinate hyperplanes 
(the minimum is over all hyperplanes and all j). 

Denote by B^~ x ,„ (- w max ) the ball of radius d min /K maiX centred at -w max in the N-1- 

dimensional hyperplane parallel to V compl , and by C^~ l j K (— w max ) the cone comprised of 
rays originating at and intersecting with the ball B^~^ K (— w max ). By linearity of maps 
Mj (and hence their superpositions) and by construction of the set C^~ l j K ^{— w max ), if 
y e C j~ 1 ik (-w max ), then 

(M^fMj-x+y G R N for all 1 < j < m and k > 

and 

lim (M u) ) k M^ 1A y = -oo for all l<j<m. 

k— >oo 

By the same arguments as employed in the proof of Theorem HI we can show that for 
any R < there exists S < such that y G C ( ^ 1 T 1 /^ max (-w max ) n U s implies 

(Af W) ) fc Af i _i, 1 y G U R for all 1 < j < m and k > 0. 
Since the measure of the set C^~^ Km ^(— w max ) n Us is positive, part (b) is proved. QED 
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5 Bifurcations of heteroclinic cycles 



In this section we consider co dimension-one bifurcations of type Z heteroclinic cycles. We 
assume that the system p]) depends on a scalar parameter, a, i.e. it has the form 

x = /(x,a), /:K"xR4l" (31) 

A bifurcation takes place at a = and the cycle exists at interval [a_,a+], where a_ < 
and a + > 0. We do not study bifurcations occurring when the cycle ceases to exist, which 
happens, e.g., if for some j a contracting, expanding or radial eigenvalue of df(£j) vanishes 
at a = 0. 

5.1 Transverse bifurcations 

Suppose a transverse eigenvalue of a steady state £j becomes positive at a = 0. In subsection 
12.41 we have proved that in a small neighbourhood of the cycle the Poincare map can be 
approximated by a superposition of maps. Near steady states only the linear part of / 
has been taken into account. The approximation is accurate in a neighbourhood of the 
steady state, provided in this neighbourhood the linear part is significantly larger than the 
omitted nonlinear terms in Taylor's expansion of /. As an eigenvalue of df(£j) tends to zero, 
the neighbourhood of £j, where the approximation is valid, shrinks. At a = 0, when the 
eigenvalue of df(£j) is zero, the collection of maps can not be used to study the behaviour 
of trajectories in the vicinity of the cycle. However, we can still comment on bifurcations 
of the cycle. Stability of the cycle can change in a bifurcation of a steady state £j, where a 
transverse eigenvalue of df(£j) becomes positive at a = 0, in the following ways: 

(a) The cycle is asymptotically stable for a < and fragmentarily asymptotically stable 
for a > 0. This happens, if conditions (i)-(iii) of Lemma H] remain satified for a > 0. 
An example of such a bifurcation is studied in the paper [6]. 

(b) The cycle is asymptotically stable for a < and unstable for a > 0. This is the case, if 
at least one of conditions (i)-(iii) is violated for a > 0. Examples of such bifurcations 
of homoclinic cycles are studied in [5] . 

(c) The cycle is fragmentarily asymptotically stable for a < and remains such for a > 
(conditions (i)-(iii) hold true for a > 0). 

(d) The cycle is fragmentarily asymptotically stable for a < and unstable for a > (at 
least one of conditions (i)-(iii) is violated for a > 0). 

When a transverse eigenvalue of df(£j) changes its sign, a pair of steady states, £j and 
emerges in a pitchfork bifurcation. The type of the global object appearing in the bifurcation 
depends on whether there exist structurally stable heteroclinic connections £j — > and 
£j ~~ ( or possibly — > with / > 1). If both connections exist, a new heteroclinic 
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cycle (£1, . . . , . . . , £ m ) is created in the bifurcation, as e.g. in the system considered 

in [B]. If such a pair of connections does not exist, more complex objects can possibly emerge, 
for instance, a depth-two heteroclinic cycle involving a connection from £J to the original 
heteroclinic cycle X. 

5.2 Resonance bifurcations 

In this subsection we assume that none of the eigenvalues of df(£j) for any j crosses the 
imaginary axis at the interval a_ < a < a + . Hence, there exists a neighbourhood of the 
heteroclinic cycle, where trajectories of the system ( )3T!) are accurately approximated by the 
superposition of maps gj, 1 < j < m. Therefore, bifurcations of the cycles can be investigated 
by studying bifurcations of fixed points of the collection of maps associated with the cycle. 
Here we consider, what happens if conditions (i)-(iii) of Lemma H] for fragmentary asymptotic 
stability of a fixed point (w, z) = of a collection of maps are broken. 

In subsection 14.21 we have divided eigenvalues and eigenvectors of a transition matrix 
into two groups, which we have called significant and insignificant. Suppose vectors in 
the basis in /T. (recall that the basis consists of eigenvectors of df(^j) ) are arranged in 
such a way that the first n s basis vectors are significant, and the last rij ones are insignificant 
(N = n s + rii). The first basis vector is the expanding eigenvector of df(^j). Significant 
eigenvalues are eigenvalues of the left upper n s x n s submatrix of M^>\ their absolute values 
are generically distinct from unity. All components of the associated eigenvectors generically 
do not vanish. The significant eigenvalues are supposed to be ordered, so that ReAj > 
ReA J+ i for 1 < j < n s . Insignificant eigenvalues are one in absolute value and the first n s 
components of the associated eigenvectors are zero. The transition matrix restricted to 
V ms is a matrix of permutation of basis vectors. A permutation is a combination of cyclic 
permutations. We arrange the insignificant eigenvectors in such a way, that each consequent 
rii vectors are involved in the same cyclic permutation of length rii and the permutation is 
v^i-i+i yJVj-i+i w h ere _/v z - Hs+ni+ _ _ _ +ril an d n . = Hl+ _ _ .+n L . 

Conditions (i)-(ii) of Lemma H] for A max are now replaced by conditions 

(i') Ai is real and Aj > 1; 

(ii') [Ai| > |Aj| for any 2<j<N. 

Denote by £ components of a vector in the basis comprised of eigenvectors of the matrix 
M. (As above, we have denoted M := M^; our arguments are applicable to for any 
j.) The map Ai, the superposition of the maps fl23|) . defines the iterates 

77 n+1 = Mr7 n :=M77 n + c. (32) 

In these coordinates the iterates take the form 

Cj.n+i = -\?Cj,n + dj (33) 
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for real Xj, and 

Ci,n+1 = a jCj,n + 0jCj+l,n + dj, Ci+l,n+l = a jCj+l,n ~ PjCj,n + (34) 

for complex Xj = otj ± i/3j . 

Lemma 5 Consider the following dynamical systems in R ((a)- (c)) and M. 2 (d): 

(a) x n+ \ = Xx n + d where X ^ ±1 is real; 

(b) x n+ i = Xx n + d, where X = — 1; 

(c) x n+ \ = Xx n + d, where X = 1; 

(d) x l n+l = ax\ + f3x 2 n + d\ x 2 n+l = ax 2 n - $x\ + d 2 . 
Fixed points and periodic orbits of these systems are: 

(a) x = d/(l — X) and x = ±oo (if X < 0, x = ±oo is a period-two orbit). The fixed point 
x = d/(l — A) is stable for \X\ < 1, and the fixed points x = ±oo are stable for \X\ > 1; 

(b) x = d/2 is a fixed point, any real number is a period-two orbit; 

(c) x = ±oo are the only fixed points. The fixed point x = oo is stable for d > 0, and the 
fixed point x = — oo is stable for d < 0; 

(e) x 1 = (d l (l -a) + f3d 2 )/((a - l) 2 + (3 2 )), x 2 = (d 2 (l -a)- (3d 1 )/ ({a - l) 2 + f3 2 )) is a 
unique fixed point, stable for \X\ < 1. 

The statements of the lemma are directly verified by a simple algebra. 

Theorem 6 Let Mj, 1 < j < N, be basic transition matrices of the collection of maps {g™}, 
associated with a type Z heteroclinic cycle. (This implies that the transition matrices have 
the form (f^P-j Suppose that 



(i) the entries bjj of the basic transition matrices depend continuously on a; 

(ii) for a_ < a < a + condition (Hi) of Lemma^ is satisfied for all ; 

(Hi) for a- < a < significant eigenvalues of matrix M satisfy the conditions Xi > 1, and 
\Xj\ < 1 for 2 < j < n s . 

(iv) for < a < a + all significant eigenvalues of matrix M satisfy \Xj\ < 1 for all 1 < j < 
n s . 

Then 

(a) For < a < the fixed point of the collection of maps {g™} is fragmentarily 
asymptotically stable, and for a > it is unstable; 
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(b) Suppose d\ < in and the point rj 1 is defined by the relations 



^ 1 

l=i;Xiis real 



X- 3 



i,l . 



E 

i=l ; A;is complex 



d l (l-Re(Xi))±d l±1 Im(Xi 
|Ai-l| 2 



-w A 



for l<j< n s ; 



(35) 



r/j = — oo for n s + 1 < j < N. 



(In the second sum the sign ± coincides with the sign of the imaginary part o/A/. The 
condition that a component is — oo implies that in the original coordinates (w, z) the 
respective component vanishes.) Then rj 1 is an asymptotically stable fixed point of the 
collection {g™} for < a < a + . Components nj (E3J) for 1 < j < n s satisfy 

lim ri] = — oo. 

J 



(c) If d\ > 0, then for a _ < a < the point |33J) is an unstable fixed point of the collection 
{g™}- Components 77] [SB) for 1 < j < n s satisfy 

lim rr) = —00. 



Proof: (a) follows from Lemma HI 

(b) Recall that insignificant eigenvectors are arranged in such a way, that each consequent 
n\ vectors are involved in the same cyclic permutation. The subspace spanned by these ni 
vectors, {v 1,7Vi+1 , v 1,iVi+2 , . . . , v 1,Ari+1 }, is an invariant subspace of map M. The action of the 
map on this subspace is presented by the matrix 



[0 


1 


. 


• °\ 








1 . 


. 








. 


. 1 


V 1 





. 


• 0/ 



(36) 



Thus, ni eigenvalues of the restriction of M into this invariant subspace are roots of unity. 
One of the eigenvalues is unity, all ni components of the associated eigenvector are equal. 

Consider the map Ai in the basis comprised of eigenvectors of M. We apply Lemma in 
each eigenspace separately to find a fixed point of the map. For eigenvalues, different from 
unity, we choose finite components in the respective eigenspaces. For the eigenvalues one 
we choose the fixed points —00. The change of coordinates from £ to rj yields the desirable 
expression flHB]) for the fixed point. 

In the sum fl35]) . the factor o? 1 /(l — Ai) is asymptotically larger than other factors involved 
in other terms. In the limit a — > +0 this factor tends to infinity, while others have finite 
limits. Hence, components rf- for 1 < j < n s tend to —00 for a — > +0, as required. 
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Next we prove that the fixed point (|35|) is asymptotically stable. Suppose that the first 
n s components, (w{ , . . . , w*' 1 ), of the eigenvector w 1 ' 1 are known (they can be found by 
calculating eigenvectors of the n s x n s upper left submatrix of M). The last rii components 
of w 1 ' 1 can be found from equations 

En s lli 1,1 

,=1^+1,^,' + W N l+ 2 

Ert s lli 1,1 

q=1 a Nl+2 , q w^ + w^ i+3 



XlW 



1,1 

N t +1 



X l W N^+2 



(37) 



where the entries of the matrix M (a separate system is obtained for each invariant 

insignificant subspace of M, and it has dimension n{). 

We solve the system (|37|) as follows. Multiplying the first equation by A™' -1 , the second 
one by A™' -2 , ... , and adding up the resultant equations we obtain 

^ «; 

w n]+i = —, ~ 5^ X T~ 3 > ] a N l+j , q w 1 /. (38) 



rn _ 7 E A l i_J E a *l+3A W q 1 - 
1 7=1 a=l 



A 

We proceed similarly to find w,' 1 , N\ + 1 < j < N t+1 . 

Since Ai tends to unity for a — > — 0, the condition it;.' > for Ni < j < iVj+i — 1 in this 
limit is equivalent to 

E E ^X' 1 > °- ( 39 ) 

3=1 q=l 

Stability of the point ( 1351) to a perturbation from the significant subspace of M follows 
from Lemma [5] and the fact that all eigenvalues in the significant subspace are less than one 
in absolute value. For a close to +0, the fixed point 77 1 has the asymptotics C(a)w 1 ' 1 , where 
C(a) < is a large constant. 

Consider now a perturbation from an invariant subspace of dimension ni of the insignifi- 
cant subspace. Suppose that D > is such that 



D > -minE^+j,^' 1 

3 9=1 



For a given 5* < 0, set R = S — niC(a)D. An initial perturbation V! from this subspace 
evolves according to the following equations: 

v£ 2 = C(a) E7=i *N l+ 2, q wY + << +3 
> Nl+ni =C(a)E n q U*N l+ni , q w^ + v^\ 
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Due to our choice of R and (I39I) . if initially all components of Vi satisfy v { < R for Ni < 
j < Ni + i, then for any k > 1 the estimate v J k < S for N t < j < iVj+i holds true. Hence, the 
point (1351) is stable in the insignificant subspace. 

To prove (c), note that the fixed point presented by expression f l35|) remains a fixed point 
of the collection of maps. It bifurcates from — oo at a = 0, and remains close to — oo for 
negative a (recall that d\ > 0). The point is unstable in the direction of w 1,1 , because the 
associated eigenvalue is larger than one. QED 

Note that if a fixed point of a collection of maps, associated with a heteroclinic cycle, 
exists near 0, then there exists a periodic orbit close to the cycle. Stability (or instability) of 
the periodic orbit follows from stability (or instability, respectively) of the fixed point of the 
collection of maps. The distance of the bifurcating fixed point from is 0(exp(D/(Ai — 1))). 
If the eigenvalue depends linearly on a (the generic case), the distance from the periodic 
orbit to the heteroclinic cycle in is proportional to e~ d// ' a ', and hence for any j the time 

required for the orbit to get from to Hj OUt ^ is proportional to l/|a|. Therefore, near the 
point of bifurcation the temporal period of the bifurcating periodic orbit behaves as 1/ 1 o; | . 

In the study of stability we have ignored matrices Aj, because only the distance to 
matters for stability. Denote A^ 1 = A^ . . . Af. Since the matrices are diagonal, the 
product is independent of the order of multiplication of the matrices. Suppose a bifurcating 
periodic orbit intersects at point x. The next intersection is at the point A ± x, coinciding 
with x if the matrix A ± is the identity, and is a different point otherwise. Therefore, the 
bifurcating orbit can have the same length as the heteroclinic cycle under consideration, or 
a twice larger length, depending on the signs of diagonal entries of matrix A^. 

In the next theorem we consider bifurcations occurring when |Ai| > 1, but some con- 
dition^) of Lemma H] fail at the point of bifurcation. In this case the fixed point of the 
collection of maps {g™} ceases to be fragmentarily asymptotically stable at the point of 
bifurcation and no new objects (fixed points or invariant sets of other types) bifurcate from 
it. Thus, the respective heteroclinic cycle looses fragmentary asymptotic stability without 
emergence of bifurcating objects. 

If entries of a matrix depend on a single parameter, generically only two real eigenvalues 
of the matrix can become equal as the parameter is varied, and at the critical parameter value 
the restriction of the matrix onto the eigenspace, associated with this eigenvalue, is a Jordan 
cell. On further variation of the parameter the two eigenvalues in this two-dimensional 
subspace change to a pair of complex conjugate ones. 

Theorem 7 Let Mj, 1 < j < m, be basic transition matrices of the collection of maps {g™} 
associated with a type Z heteroclinic cycle. (In particular, Mj have the form (24\).) Suppose 

(i) the entries bjj of the transition matrices depend continuously on a; 

(ii) for a- < a < a + the inequalities |Ai| > 1 and \j ^ 1 for 1 < j < n s are satisfied; 

(Hi) for a- < a < all conditions of Lemma^on eigenvectors and eigenvalues of matrices 
are satisfied; 



26 



(iv) for < a < a + some conditions of Lemma\4\ are not satisfied, i.e. at least one of the 
following possibilities holds true: 

(iv.l) Ai is complex, X\ = A 2 and lim^+o Im{X\) = lim Q __> +0 7m(A 2 ) = 0; 

(iv.2) there exists Xj < 0, such that \Xj\ > X± and limQ._5._1_0 |Aj| = Ai; 

(iv.S) two eigenvalues Xj and A J+ i = Xj are complex, \Xj\ > Ai and linia^+o \Xj\ = Ai; 

(iv.4) there exists q such that w^ 1 < 0, wj' 1 > for I ^ q and lim^+o w^' 1 = 0. 

Then 

(a) For o-_ < a < the fixed point of the collection of maps {g™} is fragmentarily 
asymptotically stable, and for < a < a + it is unstable. 

(b) There exists an S < 0, such that the map Ai := Ai^ does not have fixed points in Us 
( other than -00 ) for any a_ < a < a + . 

(c) There exists an S < 0, such that for any y _ Us one of the following statements is 
satisfied: 

- there exist K > and h > 1, such that |M fc+1 y| > h\M k y\ for any k > K 
or 

- for any K > there exists k > K such that M k y Us- 

(In case (c) the map M has no periodic orbits or other invariant sets in Us, other than 
-00.) 

Proof: (a) follows from Lemma HI 

To prove (b), note that by Lemma [5] the first n s coordinates of any fixed point 77 1 of the 
map Ai are -00, or otherwise they are given by ( |35l) . For S < satisfying 

S< min r,){a), (41) 

l<j<n s , __ <Q<o+ 

the map M has no fixed points in Us other than -00. (The values of 1 < j < n s , are 
bounded for a__ < a < a+, because no significant Xj is equal to unity for such a.) 

To prove (c), set S as above. Define fj 1 by relations 



Edi 1, -A d t (l - Re(Xi)) ± d l±1 lm(Xi) u 
T^x^ + ^ \x^W W ' J 

/=i;A ; is real z=i;A ; is complex 



for 1 < j < n s , 
fj] = for n s + 1 < j < N. 



(42) 
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Let y G Us he represented in the basis, comprised of eigenvectors of M, with the origin 
shifted to fj 1 . First, we note that at least one of the first n s coefficients Q in the sum 

N 
9=1 

does not vanish (otherwise by f l4Tj) the point y would be outside Us)- Let A max be the 
largest in absolute value significant eigenvalue associated with a non- vanishing coefficient. 
For k —> oo, the iterates M k y become aligned with the associated eigenvector w max (or with 
the respective 2-dimensional eigenspace, if A max is complex). If A max is real and A max > 1, 
then for large k the iterates M k y behave as (A max ) fe w max , implying that they either tend to 
-oo or escape from Us, depending on the signs of components of w max . If |A max | < 1, then 
the iterates M k y are attracted by which is outside Us due to our choice of S. If A max 
is complex, or if it is real and A max < — 1, then for any K > there exists M k y with k > K 
outside R*. QED 

The bifurcation considered in Theorem [7J was studied by Postlethwaite [17] for a particular 
dynamical system in M 4 . She proved that "there are no dynamical structures which merge 
with the cycle at the point of stability loss" , in agreement with our theorem. 



6 Homoclinic cycles 
6.1 Transition matrices 

A homoclinic cycle is a heteroclinic cycle, where all equilibria are related by a symmetry 7, 
7£j = The transition matrix of a homoclinic cycle is a basic transition matrix 

which is a product of a permutation matrix A and a local matrix B. Any permutation is 
a composition of cyclic permutations. Suppose the permutation, defined by matrix A, is a 
combination of L + 1 cyclic permutations. The first permutation involves n s significant basis 
vectors, and the last rii basis vectors are insignificant. Basis vectors can be ordered in such 
a way, that 

= ( e, + x, j * N> for any 1 < I < L 

where as above we have denoted iVj = n s + n\ + n 2 + . . . + n\. Without any restriction of 
genericity we can assume that e x is the expanding eigenvector of df(^). 

For the basis vectors ordered according to the permutation matrix A is comprised 
of L non- vanishing diagonal blocks n\ x ni, each being of the form (136|) . Therefore, the 
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transition matrix is 





I 


b 2 


1 


... 








... 


... \ 








o 


1 ... 


o 


o 


... 




















... 








bi 





... 








... 




M = 




b ni +2 





... 





1 


... 






b ni +3 





... 








1 ... 








b ni +i 





... 


1 





... 






v 




... j 



(44) 



6.2 Two examples of type Z homoclinic cycles 

In this subsection we present two general examples of type Z homoclinic cycles. In both 
cases, a foundation for our constructions is a known homoclinic cycle in W 1 with an empty 
insignificant subspace (i.e., all eigenvalues are significant and the permutation matrix A is 
just a single cyclic permutation involving n s vectors). We add n\ insignificant transverse 
directions, where admittable values of ri\ depend on n. This step (enlargement of insignificant 
subspace by adding new dimensions) can be repeated any number of times. 

The first example. Postlethwaite and Dawes [18] presented an example of system with 
a Z n k symmetry group possessing a type Z homoclinic cycle. The subg roup ^2 acts on 
¥L n as n reflections 

Sj {x\ , . . . , Xj , . . . , X n ) {x\ , . . . , Xj , . . . , x n ) , 

the subgroup Z n is generated by the cyclic permutation 

PyX\, X2-, ■ ■ ■ , 3^n) (^nj • - - > ^n—l)- 

Suppose, 

- the system ([T]) has n equilibria £j, related by the symmetry p, p£j = each equilibrium 
has an isotropy subgroup Z!?"" 1 " 2 ; 

- the unstable manifold of £j is one- dimensional, has the isotropy subgroup Z^ 1 " 1 " 1 and is 
attracted by 

Under these assumptions the system has a heteroclinic cycle. The transition matrix of the 
cycle is comprised of a block of size (n t + 1) x (n t + 1) (recall that n = n t + n r + 2). Note that 
instead of the condition that the unstable manifold of £j is one-dimensional we can impose 
the weaker condition that the connection Kj = W u (£j) PI W s (C,j + i) is one-dimensional and it 
belongs to Fix(Z2 r+1 ). 

We can enlarge the dimension of ([T|) by any K dividing n, and do it several times, in the 
following way. Consider the action of the group Z n+K x Z2 +i<: , where the subgroup Z^" 1 "^ 
acts on R n+is: as n + K reflections 

Sj(X\^ • • • , Xj) • • • , 2-71+A') (^1) • • • ) "EjJ • • • 5 ^Tl+K^l 
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the subgroup Z n is generated by the combination of cyclic permutations 

p{X\, X 2j • • • i X n , x n+2i ■ ■ ■ > •^ra-f.K') (•^n; X\, ■ ■ ■ j ■En—li X n +Ki X n +li • • • j 2'n+_ft"— l)- (45) 

As a result of this modification, the size of the transition matrix increases to (n t + 1 + if) x 
(rii + 1 + if) and a diagonal block ( l36i) of size if x if is added to the matrix A. 

A general third order system in R n+X with the above symmetry group is as follows: 

n+K 

±x = xx{v x + 2j cixf), 
1=1 

K n/K-1 K 

X n+1 = X n+1 (u 2 + ^ dl ^ tf+s + d n+l X l+l)- 
1=1 s=0 1=1 

(Equations for other Xj are obtained by application of the symmetry p.) If V\C\ < 0, the 
system has n equilibria on coordinate axes, related by the symmetry p. Sufficient conditions 
on vi and q, 1 < I < n, for existence in R™ of a homoclinic cycle connecting the equilibria 
are presented in [18]. Our construction thus implies existence of a homoclinic cycle in the 
extended (n + if )-dimensional space for all values of i/ 2 , cj, n + 1 < I < n + if anf d\. 

The second example follows [6] . The authors consider a heteroclinic cycle in R 3 [H [19] 
in a system ([T]) with the symmetry group Z 2 x Z 2 , generated by the symmetries 

Si(Xi,X 2 ,X 3 ) = (xi,x 2 , -x 3 ), 

S 2 (xi,X2,X 3 ) = (xi, ~X 2 ,X 3 ), 

p(x 1 ,x 2 ,x 3 ) = (-x 1 ,x 3 ,x 2 ). 

In [6], the dimension of the system is increased up to 5 and the symmetry group is 
enlarged to Z 4 x Z 3 , by adding the symmetry 

s 3 (xx,x 2 , x 3 , x 4 , x 5 ) = (xx,x 2 , x 3 , x 4 , -x s ). (46) 

Symmetries Si and s 2 act on the added dimensions {x±,x$) trivially, and p is modified to 
become 

p(xi,x 2 ,x 3 ,X4,x 5 ) = (-x 1 ,x 3 ,x 2 ,-x 5 ,x i ). 
The system in R 5 with the above symmetry group truncated at the third order is as follows: 

±1 = V X X\ + ci(x\ - x\) + c 2 (x\ - x\) + xi(c 3 x\ + c 4 (x 2 + x\) + c 5 (x| + x\)) 

X 2 = X 2 {v 2 + C 6 Xi + C 7 x\ + C$xl + C 9 xl + C\qx\ + C\\xf) 

x 3 = x 3 {y 2 - c 6 x 1 + c 7 x\ + c 8 x 2 + C9X 2 + C10X 2 + cux\) (47) 

X4 = x 4 (z/ 3 + diXi + d 2 x\ + <i 3 X2 + d±x\ + d^x\ + <i 6 x 2 ) 

X4 = x 5 (z/ 3 — d\X\ + d 2 x\ + d 3 x\ + d±x\ + <i 5 a; 2 + d 6 x 2 ). 
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If < 0, the system posesses two equilibria on the x axis related by the symmetry p. 
Sufficient conditions for existence of a homoclinic cycle in M. 3 involving the equilibria are 
presented in [H US]- In [S] this homoclinic cycle in M 5 is studied for particular values of 
coefficients u, c and d in (jlTf) . 

By the same procedure dimension of ([I]) can be incremented by 1 or 2 any number of 
times. Also, the dimension the system can be expanded by choosing dimension of Lj larger 
than unity, like in the first example considered in this subsection (the resultant transition 
matrix is not modified). 

6.3 Stability of a cycle when all transverse eigenvalues are nega- 
tive 

If all transverse eigenvalues of df(£) are negative, the cycle can be unstable or asymptoti- 
cally stable (see theorem 0j). Suppose vectors comprising the basis are ordered according to 
(j43p . The eigenvalues associated with eigenvectors from the insignificant subspace are one 
in absolute value (see subsection 14. 2p . As proved in [TS], the dominant eigenvalue of the left 
upper n s x n s submatrix is real and larger than one, if and only if 

61 + 62 + ... + ^. >1. (48) 

By Theorem [TJ condition (1481) is necessary and sufficient for a homoclinic cycle of type Z, 
whose all transverse eigenvalues are negative, to be asymptotically stable. 

6.4 Two simple cases 

If the dimension of the significant subspace is one or two, then the dominant eigenvalue and 
the associated eigenvector can be explicitly calculated in the terms of eigenvalues of df(£). 
The left upper submatrix is then 



(£»i) or 




for the dimension one or two, respectively. 
If the dimension is one, the eigenvalue is 

A = h, 

and the necessary condition for stability is hi > 1. If the dimension is two, the dominant 
eigenvalue is 

b 2 + ^/b\ - 4&! 
A " 2 ' 

the necessary conditions for stability are 62 > and b\ + 62 > 1 (see p3]). The conditions 
imply, that the first two components of the associated eigenvector have the same sign. We 
assume that they are positive. 
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Let us calculate other components of the eigenvector w of the transition matrix, associ- 
ated with the dominant eigenvalue A. (For simplicity, upper indices are omitted.) Consider 
components (wfy^+i, iyjv ! _ 1 +2, • • • , Wjv,) involved in a single permutation cycle of length n\. 
The components satisfy the equations 

b Nl _ 1+ iWi + w Nl _ 1+2 = \w Nl _ 1+ i, 

&AT ! _ 1 +2^1 + W N{ _ 1+3 = Xw Nl _ 1+2 , 

b Nl w 1 + w Nl _ 1+1 = Xw Nl . 

Multiplying the equations by X h ( k+1 ' ni \ . .. ; \ h ( k + n i> n i) ^ where h(K,ni) = mod ni (K), re- 
spectively, for < k < ni — 1, and adding them up we obtain 

{b Nl _ 1+1 X h ^ k+1 ^ + ... + b Nl X h ^ n ^) Wl = (A n « - l)w Nl _ 1+h . 

Thus, the components Wjvj_i+1; • • • , w Ni are of the same sign as Wi as long as 

b Nl _ 1+1 X h(k+1 ' ni) + ... + b Nl X h{k+nuni) > for any < k < m - 1. (49) 

If the length of the first permutation cycle is larger than two, the same conditions guar- 
antee stability of the cycle: conditions (i)-(iii) of Lemma H] for eigenvalues and eigenvectors 
of the significant n s x n s submatrix, and condition (TJ9]) for each insignificant permutation 
cycle. However, we can not express the eigenvalue A in the terms of eigenvalues of df(£), 
except for n s = 3 and 4, but in these cases the expressions are too complex to be of practical 
use. 

7 Conclusion 

We have defined type Z heteroclinic cycles as simple cycles with a certain action of subgroups 
of the symmetry group of the dynamical system. We have introduced the notion of frag- 
mentary asymptotic stability of an invariant set and have derived necessary and sufficient 
conditions for asymptotic stability or fragmentary asymptotic stability of type Z heteroclinic 
cycles. For a type Z cycle we have calculated basic transition matrices, which depend on 
eigenvalues of linearisations near steady states comprising the cycle and the action of the 
system symmetry group. 

We have proved that a type Z cycle is asymptotically stable, if all transverse eigenvalues of 
linearisations near the steady states are negative and the largest in absolute value eigenvalue 
of transition matrices (which are products of the basic transition matrices; all transition 
matrices are similar) are larger than one in absolute value. For such type Z cycles, eigenvalues 
of transition matrices play the same role, as eigenvalues of the linearisation in the study of 
stability of steady states and eigenvalues of the linearisation of the Poincare map in the 
study of stability of periodic orbits. 
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Type Z cycle, where some of transverse eigenvalues are positive, can be fragment arily 
asymptotically stable. We have derived a criterion for fragmentary asymptotic stability in 
the terms of eigenvalues and eigenvectors of its transition matrices. 

We have studied bifurcations occurring when conditions for fragmentary asymptotic sta- 
bility cease to be satisfied. Two types of bifurcations have been identified: transverse and 
resonant ones. A detailed study of transverse bifurcations is problematic, because the col- 
lection of maps, that we use to study stability, does not approximate correctly trajectories 
of the system for the critical parameter value. However, we comment how stability of the 
cycle can change in a bifurcation of a steady state, where a transverse eigenvalue of the lin- 
earisation near a steady state vanishes. For a resonant bifurcation, we prove that either it is 
accompanied by a birth of a periodic orbit (if the dominant eigenvalue of a transition matrix 
becomes smaller than one), or no new invariant objects emerge in it (if other conditions for 
fragmentary asymptotic stability are broken). 

We anticipate the following continuations of the study. 

Apparently, transition matrices can be used to study stability of simple cycles of other 
kinds - however, it is not evident how to define a transition matrix for more complex cycles. 

We have not considered all issues related to asymptotic stability and bifurcations of type Z 
heteroclinic cycles. As we have noted in sectional bifurcations resulting in destruction of the 
cycle are yet to be investigated. Our study of transverse bifurcations is by far incomplete. 
We have assumed that in a neighbourhood of the cycle the behaviour of trajectories is 
accurately approximated by a collection of maps, in which only linearisations are taken 
into account. This is true in a neighbourhood, where the linear part is significantly larger 
than the nonlinearity, but when a (transverse) eigenvalue of the linearisation vanishes, the 
neighbourhood of a steady state where this holds true shrinks to zero. For a more detailed 
analysis of the bifurcation, some non-linear terms in the local expansion of the r.h.s. / should 
also be included into the local map. 

In subsection 16.21 we have presented two general examples of type Z homoclinic cycles. 
The question arisies whether type Z homoclinic cycles of other kinds exist; if they do exist, 
their classification is wanted. (Heteroclinic cycles are apparently too diverse for any such 
classification to be useful.) 
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